Stability of solitons in PT-symmetric couplers 
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Families of analytical solutions are found for symmetric and antisymmetric solitons in the dual-core system 
with the Kerr nonlinearity and PT-balanced gain and loss. The crucial issue is stability of the solitons. 
A stability region is obtained in an analytical form, and verified by simulations, for the PT-symmetric 
solitons. For the antisymmetric ones, the stability border is found in a numerical form. Moving solitons 
of both types collide elastically. The two soliton species merge into one in the "supersymmetric"case, 
with equal coefficients of the gain, loss and inter-core coupling. These solitons feature a subexponential in- 
stability, which may be suppressed by periodic switching ("management"). © 2011 Optical Society of America 
OCIS codes: 060.5530, 190.6135, 190.5940, 230.4320 



Dissipative media featuring the parity-time (VT) sym- 
metry have recently drawn a great deal of attention. 
The introduction of this symmetry in optics followed 
works extending the canonical quantum theory to non- 
Hermitian Hamiltonians that may exhibit a real spec- 
trum [1]. The Hamiltonian is PT-symmetric if it in- 
cludes a complex potential V{x) which satisfies con- 
straint V{x) = x). Such potentials were proposed 
[2]- [7] and realized [8, 9] in optics, by juxtaposing spa- 
tially symmetric patterns of the refractive index and ap- 
propriately placed gain and loss elements, see Ref. [7] for 
a review. Nonlinear VT systems [10]- [15] and respective 
solitons [16] were introduced too. 

A medium which is akin to VT systems is a dual- 
core waveguide with gain and loss acting separately in 
two cores, which are linearly coupled by the tunneling 
of light [17]. This system predicts stable ID solitons in 
optical [17]- [19] and plasmonic [20] waveguides with the 
Kerr nonlinearity, as well as 2D dissipative solitons and 
vortices [21,22]. The system is made VT symmetric by 
adopting equal strengths of the gain and loss in the cores. 
A challenging problem is the stability of solitons, as sta- 
ble pulses in the dual-core system were previously found 
far from the point of the VT symmetry [18-20]. We pro- 
duce two families of exact soliton solutions for the VT- 
symmetric system, which correspond to symmetric and 
antisymmetric solitons in the ordinary dual-core cou- 
pler [23]- [26]. For the former family, an exact stability 
border is found analytically, and verified by simulations. 
For the PT-antisymmetric solitons, the stability region 
is identified in a numerical form. In the "super symmet- 
ric" limit, when the gain and loss coefficients coincide 
with the constant of the inter-core coupling, both fami- 
lies merge into solitons which are subject to a subexpo- 
nential instability. We demonstrate that this instability 
can be suppressed by means of a "management" tech- 
nique, i.e., periodic switching of the gain, loss, and cou- 
pling. 

The transmission of light or plasmons in the dual-core 



waveguide is described by the linearly coupled equations 
for amplitudes u{z,t) and v{z,t) in the active and passive 
cores [17]- [20]: 

iuz -\- {l/2)utt -\- \u\'^u — iju -\- K,v = 0, 
iVz^{l/2)vu^\v\'^v^iTv^Ku = 0, (1) 

where z is the propagation distance and t the reduced 
time or transverse coordinate in the temporal- or spatial- 
domain system. Coefficients accounting for the disper- 
sion or diffraction, Kerr nonlinearity, and inter-core cou- 
pling, are normalized to be 1, while 7 and F are coef- 
ficients of the linear gain and loss in the two cores. The 
VT symmetry in Eqs. (1) corresponds to F = 7. Strictly 
speaking, the PT-balanced gain and loss correspond to 
the border between stable and unstable systems: The 
zero solution, u = v = 0, is unstable at 7 > F, while the 
stability region was found at 7 < F [17, 18]. 

Obviously, any solution to the nonlinear Schrodinger 
(NLS) equation (with a frequency shift), iUz-\-{l/2)Uu-\- 
\U\'^U ± ^/l — 7^[/ = 0, gives rise to two exact solutions 
of the PT-symmetric system, provided that 7 < 1: 

v= (^i-f± \/r^^) u = U{z, t) (2) 

(recall we fix = 1). For 7 = 0, solutions (2) with + and 
— amount, respectively, to the symmetric and antisym- 
metric modes in the dual-core coupler [23]- [26], therefore 
we call the respective solutions (2) PT-symmetric and 
PT- antisymmetric ones. In the limit of 7 = 1, exact so- 
lutions (2) reduce to a single one, v = iu = U{z,t). In 
particular, PT- symmetric and antisymmetric solitons, 
with arbitrary amplitude 77, are generated by the NLS 
solitons, 

U {z, t) = r] exp ± ^/l - 7^^ z^ sech {r]t) . (3) 

As concerns stability of the solitons, it is again relevant 
to compare the T^T-symmetric system to its counterpart 
with 7 = 0. The exact result for the latter system is that 
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the symmetric solitons are unstable against symmetry- 
breaking perturbations at rf > r]'^^,^ (7 = 0) =4/3 [23]. 
The antisymmetric sohtons are unstable against oscilla- 
tory perturbations at r]'^ > fj^g^^ ~ 0.8 [25]. The latter 
instability threshold was found in a numerical form. 

The analysis of the instability against an- 
tisymmetric perturbations, du = —dv, which 
leads to the exact result rj'^s,^ (7 = 0) = 4/3, 
can be extended for the T^T-symmetric solitons. 
The respective perturbation 5u at the critical 
point, T]'^ = 7?max5 obcys the linearized equation, 
{4^/1 - 7' - d^dt^ + [1 - 6sech' (r/maxt)] }5u = 
0. This solvable equation yields 

^Lx(7) = (4/3)x/r^. (4) 

We have performed simulations of the evolution of per- 
turbed solitons within the framework of Eqs. (1), aiming 
to verify the stability border (4) of the T^T-symmetric 
solitons, and identify a border for their antisymmetric 
counterparts. Perturbations were introduced by indepen- 
dently changing the initial amplitudes of both compo- 
nents by ±3%. The results are summarized in Fig. 1 
for both soliton families. The numerical stability bor- 
der for the symmetric solitons goes somewhat below 
the analytical one (4) because the finite perturbations 
used in the simulations are actually rather strong. Ac- 
cordingly, taking smaller perturbations, we obtain the 
numerical stability border approaching the analytical 
limit. For instance, at 7 = 0.5, amplitude perturbations 
±5%, ±3%, and ±1% give rise to the stability border at 
^max = 1-02, 1.055, and 1.08, respectively, while Eq. (4) 
yields r]'^^^ ~ 1.15 in the same case. 
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instability of solitons with r] > rjmax leads to a blowup 
of field u and decay of v in the direct simulations (not 
shown here). 




Fig. 2. (Color online) Examples of the relaxation of sta- 
ble perturbed T^T-symmetric (a) and antisymmetric (b) 
solitons, with 7 = 0.5,/^ = 0.9 and 7 = 0.5,/^ = 0.7, 
respectively. Here and in the case shown in Fig. 3, the 
evolution of the ^'-component is similar to what is shown 
for u{z^ t). 

The Galilean invariance of Eqs. (1) suggests to con- 
sider collisions between stable solitons, setting them in 
motion by means of boosting^ {u^^} {u,v}ex.p{±ixt), 
with frequency shift x- The results is that the colli- 
sions always seem elastic. A typical collision between 
PT-symmetric and antisymmetric solitons is displayed 
in Fig. 3. 




Fig. 1. (Color online) The analytical stability border 
(4) and its numerically found counterpart for the VT- 
symmetric solitons, and the numerical stability border 
for the 7^ T- antisymmetric solitons. The solitons with 
amplitude r] are stable at 7^ < 77max- 

Examples of the stable propagation of perturbed soli- 
tons of both types are displayed in Fig. 2. These ex- 
amples are shown for the solitons located close to the 
stability border, cf. Fig. 1 [in particular, Eq. (4) yields 
7/cr(7 = 0.5) ^ 1.07, cf. T] = 0.9 in Fig. 2(a)]. Unlike 
the dual-core coupler [17]- [26], the VT system cannot 
support asymmetric solitons, as the balance between the 
gain and loss is impossible for them. Accordingly, the 
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Fig. 3. (Color online) The elastic collision between stable 
solitons of the two different types, with 77 = 0.7, boosted 
by frequency shift x = at 7 = 0.7. 

In the case of the "supersymmetry" , 7 = 1, the sta- 
bility regions in Fig. 1 shrink to nil. In this case, the 
linearization of Eqs. (1) with 7 = F = around NLS 
solution (2) leads to the following equations for pertur- 
bations Su and Sv: 

L {Su + iSv) =0,1/ {Su — iSv) = —2m {Su + iSv) , (5) 

where LSu = [id, + {l/2)dtt + 2\U\^] Su^U^Su"" defines 
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the linearization operator of the NLS equation. If the 
NLS solution is stable, the first equation (5) produces no 
instability, while the second gives rise to the resonance^ 
as {Su + i5v) is a zero mode of entire operator L, which 
includes term idz . According to the linear-resonance the- 
ory, the respective perturbation {Su — iSv) will be un- 
stable, growing ~ z (rather than exponentially). Direct 
simulations of Eqs. (1) with 7 = confirm that the soli- 
tons are unstable, the character of the instability being 
consistent with its subexponential character. 

The "supersymmetric" solitons may be stabilized by 
means of the management technique [27,28], periodically 
reversing the common sign of 7 = F = which does 
not break the supersymmetry of system (1). Flipping 7 
and F means switch of the gain between the two cores, 
which is possible if both are doped [29]. The coupling 
coefficient, hz^ cannot flip by itself, but the signal in one 
core may pass 7r-shifting plates, which is tantamount to 
the periodic sign reversal of k,. 

Ansatz (2) still yields exact solutions to Eqs. (1) with 
coefficients F = 7 = subjected to the periodic manage- 
ment. On the other hand, the replacement of K: by the 
periodically flipping coefficient destroys the resonance in 
Eq. (5). In simulations, this management mode indeed 
provides for self-trapping of robust solitons from inputs 
significantly different from the exact soliton solution. An 
example is shown in Fig. 4 for the input with a relatively 
large perturbation, in the form of a 10% amplitude devi- 
ation from the exact solution. A detailed analysis of the 
management will be reported elsewhere. 




Fig. 4. (Color online) Self-trapping of the soliton in the 
"supersymmetric" system (F = 7 = = 1) subject to 
the management with period Az = 2. The input is uq = 
—ivo = 0.9sech(t). 

In conclusion, we have found the families of exact sym- 
metric and antisymmetric solitons in the T^T-symmetric 
dual-core system with the Kerr nonlinearity. For the 
symmetric family, the stability border (4) was found in 
the analytical form, and corroborated by simulations. 
The stability region for the antisymmetric solitons was 
obtained in the numerical form. Collisions between soli- 
tons of both types are elastic. In the "supersymmetric" 
limit, with the equal gain, loss and inter-core-coupling 
coefficients, the two soliton species merge into one, which 
is subject to the subexponential instability, and can be 
stabilized by means of the periodic management. It may 



be interesting to extend the spatial-domain model by in- 
cluding a periodic linear or nonlinear potential, which 
may give rise to other soliton modes [30] , and to include 
the intermode dispersion in the temporal domain [31]. 
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